Bounds on Information Propagation in Disordered Quantum Spin Chains 
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We investigate the propagation of information through the disordered XY model. We find that all 
correlations, both classical and quantum, are exponentially suppressed outside of an effective lightcone 
whose radius grows at most logarithmically with |t[. 
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How fast can information propagate through a locally in- 
teracting system? For classical systems an essentially uni- 
versal answer to this question is that the velocity of infor- 
mation propagation is bounded (often only approximately) 
by an effective speed of light. It is a more subtle issue 
to formulate equivalent velocity bounds for quantum sys- 
tems because they can encode quantum information in the 
form of qubits and therefore might be able to exploit quan- 
tum interference to propagate information faster. However, 
for quantum spin networks this is not the case: the Lieb- 
Robinson bound limits the velocity at which correlations 
can propagate Iffl . 

The Lieb-Robinson bound implies that there is an effec- 
tive light cone for two-point dynamical correlations, i.e., 
apart from an exponentially suppressed tail, two-point cor- 
relations propagate no faster than a speed of light. Sim- 
plified and alternative proofs of the Lieb-Robinson bound 
have been subsequently discovered J21 EJ 01 [5|. More re- 
cently, it has been realised that the Lieb-Robinson bound 
is strong enough to bound not only the propagation of two- 
point correlations but of any local encoding of information: 
regardless of the encoding no information (either quantum 
or classical) can propagate faster than the speed of light for 
the system J6). 

There are many consequences of the Lieb-Robinson 
bound. Apart from the aforementioned bounds on the ve- 
locity of information propagation, it has been realised that 
the Lieb-Robinson bound can be used to provide a method 
to efficiently simulate the properties of low-dimensional 
spin networks HJ |U [lOj QT). Additionally, using the 
Lieb-Robinson bound, dynamical entropy area laws for 
quantum spin systems can be obtained ll6l lT2ll . 

While the Lieb-Robinson bound is extremely general — 
it relies only on the ultraviolet cutoff imposed by lattice 
structure — it is, as a consequence, relatively weak. Thus, 
it is extremely desirable to develop stronger bounds con- 
straining the propagation of quantum information through 
systems where more is known about the structure of the in- 
teractions. One setting where one would expect stronger 
bounds to be available is that of a quantum spin system 
with disordered interactions. Such systems have attracted 
a large amount of interest as they can exhibit the striking 
phenomenon of Anderson localisation lfT3l . which means 
that information is essentially frozen: a quantum particle 
placed anywhere within a localised system diffuses only 



site 



FIG. 1: Schematic illustration of the conjectured logarithmic 
light cone for disordered systems: as time progresses information 
is exponentially attenuated outside of a light cone whose radius 
grows at most logarithmically with time. 



slightly, even for extremely large times. Thus, exploit- 
ing the parallels between bounds on information propaga- 
tion and Lieb-Robinson bounds, we are motivated to con- 
jecture that interacting spin systems with disordered inter- 
actions satisfy stronger bounds on correlation propagation 
(see Fig.[T]). More specifically, we conjecture that for quan- 
tum spin networks with disordered interactions all correla- 
tions, both quantum and classical, are suppressed outside 
of a light cone whose radius grows at most logarithmically 
in time. (Contrast this with the light cone supplied by the 
Lieb-Robinson bound: it has a radius which grows linearly 
with time.) 

In this Letter we study a setting where the dynamics of a 
class of disordered interacting spin systems can be shown 
to satisfy our logarithmic light-cone conjecture; we study 
the XY spin chain with disordered interactions in a disor- 
dered magnetic field and show that information, and hence 
correlations, are exponentially attenuated outside of a light 
cone whose radius grows logarithmically with time. The 
main result of this Letter, the logarithmic light cone for the 
disordered XY model, can by summarised by the follow- 
ing bound on the dynamic two-point correlation functions: 
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where A, and B k are local operators acting nontrivially 
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only on spins j and k respectively, Z max is the localisation 
length of our system, and c and v are constants. We ap- 
ply our new bound to study the structure of the propagator 
for large times and the scaling of the entropy of a block of 
spins in the evolving system. As a consequence, we prove 
the entropy saturation numerically observed by De Chiara 
et. al. [22 1. Our results also constitute a proof of a conjec- 
ture raised in lfT4l : namely, if two parties, Alice and Bob, 
have access to a bounded region at either end of the chain, 
respectively, then it is impossible for Alice to send any in- 
formation to Bob, regardless of how Alice encodes the in- 
formation in the single- and higher-excitation sectors. 

We consider a one-dimensional chain of n spin- 1/2 par- 
ticles with XY-model type interactions between nearest- 
neighbouring spins in an additional transverse field (e.g. a 
magnetic 5-field). We allow the coefficients of the cou- 
plings and the transverse field strength to vary from site to 
site within the spin chain. Thus, we study the evolution of 
the chain under the Hamiltonian 
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where fj,j and Uj are drawn from probability distributions 
P M and IP„ respectively, and where cr? (a E {X, Y, Z}) is 
a Pauli spin operator acting on the spin at site j. Typically, 
fij = — J for all j, however this is not necessary and we 
deal with the more general case here. 

We solve this system using the Jordan-Wigner transform 
(for an introduction to the Jordan-Wigner transform see 
|[T5ll ") which, when combined with some exact results from 
the theory of localisation, allows us to bound the dynam- 
ics of our spin chain. We note that it is not immediate that 
the dynamics of the disordered XY model satisfy a loga- 
rithmic light cone: while the XY model is equivalent to a 
free fermion model which is the original Anderson model, 
and is localised, the Jordan-Wigner transform is a nonlocal 
operation and might confound the localisation occurring in 
the free fermion system. 

Let's start by applying the Jordan-Wigner transform, 
which maps a system of interacting qubits into a system 
of free fermions. The Jordan-Wigner transform defines an- 
nihilation operators a } = (erf • • • <y^_^)a^ (where a 3 = 
|0) (1| acts on site j) and the corresponding creation op- 
erators at, which satisfy the canonical fermionic anticom- 
mutation relations. Using this we are able to rewrite the 
system Hamiltonian as H = ^™ k=1 Mj k a^a k , where the 
tridiagonal matrix M is defined via M^ k = 1\x k b^ k +\ + 

It is now possible (following the method described in 
031) to diagonalise H. After doing so we find the dynam- 
ics for the annihilation operators in the Heisenberg picture, 
with a,j(t) 
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a i^) = /~2 v 3k(t)a k , 
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where Vj k (t) = (e We now concentrate on 

bounding the Vj k {t), which in turn bounds the dynamics 
of the system. 

To proceed further, we must be more specific about our 
model: we suppose that P„ is a Cauchy distribution with 
parameter (or "width") 5 

7T [Uj — V) 2 + 5 2 

We note that, as pointed out in |[T4l . the Uj can have a dif- 
ferent distribution without affecting the qualitative aspects 
of the following results (i.e. with a different distribution we 
still get localisation of eigenstates), although of course the 
quantitative results will differ. 

Under our assumption about the distribution of the Vj, 
we know from the theory of localisation lfl6l that the eigen- 
states of the matrix M are, with high probability, exponen- 
tially localised. The matrix M is an n x n matrix (it is the 
Hamiltonian for the single-particle sector). We write {| j)} 
for the basis induced by the notation M jfe = (j\M\k). We 
can now describe explicitly how localisation manifests it- 
self in the eigenstates of M. If \E a ) is an eigenstate of M 
with eigenvalue E a , then 
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where j a is the site around which \E a ) is localised, N a = 
I (ja\E a ) | is a normalisation constant, and X a = l/l a is 
the inverse of the localisation length of eigenstate \E a ): 
small l a means that \E a ) is a highly localised eigenstate 
whilst large l a means that \E a ) is a weakly localised eigen- 
state. 

Using the localisation of the eigenstates of M we 
are able to bound the matrix elements Vj k {t): we 
see that v jk (t) = (j\e- ltH \k) = (j\k(t)) = 
ELi (E a \k) e~ itE - OW.andso 
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Defining Z max = sup Q l a and using the triangle inequality 
together with N a = \ (j a \E a ) \ < 1, we find that 
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We note that for small 5, L 



where /x is the geomet- 
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ric mean of the fj,j |[T4l . In particular, if fj,j = — J for all j, 
then fj, = J. 

Recall that v ]k = (e~ lMt )j,k- The inequality Eq. ^ is 
a quantitative statement of the result that the modulus of 
the diagonal matrix elements of e~' lMt are large, while the 
modulus of the off diagonal matrix elements decay expo- 
nentially with distance from the diagonal. This means that 
dj(t) is effectively a linear combination of only a small 
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number of a k operators — namely those for which \j — k\ 
is small. It is this fact which inhibits the spread of opera- 
tors on the chain, giving rise to the logarithmic light cone 
that we derive below. 

We now turn to the proof of the improved Lieb-Robinson 
bound for our system. We begin by bipartitioning the spin 
chain into two sections, A and B, where we assume the 
boundary between partitions is between spins m and m+1. 
We then attempt to write e ltH as a product of e ltHA and 
e ltHB . Clearly this won't be exact and so we introduce 
an operator V{t) which bridges the boundary between A 
and B, and which is designed to compensate for any errors 
introduced: 

e ltH = e tt(HA+HB) V(t). (8) 

The operator V{t) acts nontrivially on all spins in the 
chain, however, we now show that V{t) can be well ap- 
proximated by another operator, which we call V it), 
which acts only on a small number of spins. The reason 
we can do this is that V (i) acts strongly on spins which are 
close to the boundary and progressively weaker on spins as 
we move away from the boundary. To prove this approxi- 
mation is valid, we use the following differential equation 
forF(i): 

^ = iV{t)h m {t), (9) 

where h m {t) = e~ ltH h m e ltH and h m is the interaction 
term in the Hamiltonian which bridges the boundary. We 
let 0, denote a set of \£l\ spins centred on the boundary 
between the partitions A and B. We also define h^ n {t) = 
e -itfln^ me ztflY2 w here Hq contains only those interactions 
in H which act on sites in £1. We then define V Q (t) via 

d 



dt 



V n {t)=iV n (t)hl{t). 



(10) 



Clearly the operator V (t) acts nontrivially only on £1. 
The error between V{t) and V it) is bounded by 

r\t\ 

\\V{t)-V n (t)\\< \\h m {s)-hl{s)\\ds. (11) 



Calculating ||/t m (t) — 7^(^)11 15 a lengthy but straightfor- 
ward task, and we begin by using the Jordan Wigner trans- 
form to write this quantity in terms of the dj operators: 
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When we calculate ||/i m (t) — /^(t)|| we'll have to deal 
with terms such as \\al n (t)a m+1 (t) - a^(t)a^ +1 (t) ||, 
which can be bounded as follows 

\\al(t)a m+1 (t) - a%(t)a° +1 (t)\\ < 

I K> (0|| \\a m +i(t) - a^ +1 (t)\\ 

+ \K(t)-a%(t)\\\\a% l+1 (t)\\. (13) 



Now ||a m (i)|| = || a 
the problem of bounding 
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1, so we've reduced 
- h2,(t)\\ to bounding 



\a m {t) — a" (*)!!• Tne operator a" (i) is given by 
<W = ^v mk {t)a k . 
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Hence we arrive at 
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where we've used our bound Eq. ^ on \vj k \ and the fact 
that ||cifc|| = 1. 

Since Q is a set centred on the boundary between parti- 
tions A and B of the chain, we have that \m — k\ > 1 0| /2 
for all k ^ f2. Hence 



\a m (t)-aZ(t)\\<n(n-\n\)e-^ 



(16) 



and so we are finally able to conclude that 

\\h m (t) - hg(t)\\ < cn 2 e-N/ 2 '— and that 



\V(t) -V n (t)\\ < c\t\n 2 e 
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for some constant c. In particular, given e > 
0, choosing |fi| > 2/ max log (c|t|re 2 /e) ensures that 
\\V(t) - V n (t)\\ < e. That is, given any e > we can 
choose £1 to be a large enough set such that V n (t) approx- 
imates V(t) to within e. This enables us to write 
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e it{H A+ H B )yU^ + ^ 



(18) 



Following [7 1 we recursively apply the above partitioning 
procedure to find e ltH = Q(t) + 0(e), where 



Q(t) 




3=1 




(19) 



and where V\ = {Qj} is a partition of the chain into ^ 

blocks, each containing |fi| spins and where Vi = {Q' k } 
is a partition of the chain obtained by shifting V\ along by 
sites (note that £l' Q and £l' n /\ n \ are half-size blocks of 
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4^ sites each). This is our fundamental structure result for 
the dynamics of the disordered XY spin chain. 

A Lieb-Robinson bound is an upper bound on quanti- 
ties such as || [A, B(t)] ||. We now show how the above 
structure result implies a version of the Lieb-Robinson 
bound which is substantially stronger than the original. De- 
fine B{t) to be the operator which arises when we evolve 
B according to the approximation Q(t) of e ltH , namely, 

B{t) = Qit)BQpit). JThis enables us to write Bit) = 
B(t)+0(e). Note that B(t) acts trivially on all sites which 
are a distance greater than 3|0|/2 away from those sites on 
which B acts. Thus, if d(A, B) > 3|0|/2, where d(A, B) 
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is the distance between A and B, then [A, B(t)\ = 0, and 
so for a given 
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(21) 



where c and k are constants. This is the logarithmic light 
cone for the two-point dynamical correlation functions. 
Compare this to the original Lieb-Robinson bound, which 
reads 



\\[A,B(t)}\ 



< ce ki\t\ e -k 2 d(A,B) 



(22) 



To conclude we'd like to mention two consequences of 
the logarithmic light cone for the disordered XY model. 
The first is a proof of the conjecture that two parties, Al- 
ice and Bob, with access to only bounded regions A and B 
at either end of the chain, respectively, cannot use the dy- 
namics of the disordered model to send information from 
Alice to Bob. We follow the argument of 16), appropriately 
modified to take account of our stronger bound. 

Let C = L \ (A U B), where L is the chain, be the 
region that Alice and Bob cannot access. The most general 
way Alice can encode her message is via a set of unitary 
operators {U\ \k = 1,2, ... , m} on her system, where 
k is varied according to the message she wants to send. 
After a time t has elapsed the system has evolved from an 



initial state p Q to p{t) = e 



-iHt 
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We interpret this 



as a quantum channel with input p k ABC = U\pqU\ and 

output p k B {t) = tr Ac (UX(t)p UX\t)). As argued in 0, 
the output states are all very close together, as measured in 
trace norm: 

\\p B (t)-p B mi<cn 2 \t\e^ A - B ^-, 

where p B (t) = tT AC (e~ iHt p e iHt ). 

If Alice applies the unitaries {UX} according to the 
probability distribution {pk}, the amount of information 
that is sent through the channel is given by the Holevo ca- 
pacity: 



x(t) 




E^ 5 '(^(i))> 



k = l 



where S(-) is the von Neumann entropy. Applying Fannes 
inequality ifTTI we find that 

X (t) <2e(|B|-log 2 (e)), 

where e = cn 2 \t\e~ vd ^ A ' B ^ lmax . That is, Bob has to wait 
an exponentially long time (in d(A, B)) before a nontrivial 
amount of information can arrive. The optimal encoding 
for Alice to adopt was investigated in |[T8l and lfT9ll and 
completely solved in the single-use case in EDI . 

The second consequence of the logarithmic light cone 
bound is that the entropy of any contiguous block B of 



spins in a dynamically evolving product state \ip{t)) = 
e ltH 1 00 • • • 0) is bounded. Indeed, applying the argument 
of 021ED, we find that S(p B (t)) < ci + c 2 log 2 (n|t|) as 
\B\ — > oo, where Ci and c 2 are constants. This provides 
a theoretical explanation for the phenomenon numerically 
observed by De Chiara et. al. Il22l . It seems nontrivial to 
adapt the argument of |6) to prove the same result because 
their proof can't be simply modified to make use of the 
presence of disorder. 
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